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Adiabatic Faraday effect in a two-level Hamiltonian formalism 
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The helicity of a photon traversing a magnetized plasma can flip when the B-field along the 
trajectory slowly reverses. Broderick and Blandford have recently shown that this intriguing effect 
can profoundly change the usual Faraday effect for radio waves. We study this phenomenon in 
a formalism analogous to neutrino flavor oscillations: the evolution is governed by a Schrodinger 
equation for a two-level system consisting of the two photon helicities. Our treatment allows for a 
transparent physical understanding of this system and its dynamics. In particular, it allows us to 
investigate the nature of transitions at intermediate adiabaticities. 

PACS numbers: 41.20.Jb, 42.25.Lc, 95.85.Sz, 14.60.Pq 



I. INTRODUCTION 

Different polarization states of electromagnetic waves 
propagating in media and/or external -B-fields usually 
have different refractive indices, leading to a nontrivial 
evolution of a given polarization state. In astronomy, 
the rotation of the plane of polarization caused by a 
magnetized medium along the line of sight (Faraday 
effect) is the most important example. Such phenomena 
are perfectly analogous to particle oscillations and, in 
particular, to neutrino flavor oscillations where the role 
of polarization is played by flavor. 1 Moreover, in 
addition to mixing photons of different polarization, 
-B-fields also mix photons with certain other particles, 
notably gravitons, neutral pions, or hypothetical axion- 
like particles, which in this context play the role of 
additional photon polarization states 

If the medium (taken to include external fields) varies 
along the photon trajectory, the refractive indices of two 
polarization states may cross over. In the adiabatic 
limit a true crossing is avoided and opposite helicity 
states are adiabatically connected, leading to complete 
transformation. In neutrino physics, this phenomenon 
is known as the Mikheev-Smirnov-Wolfenstein (MSW) 
effect, allowing for large flavor transformations even when 
the vacuum mixing angle is small Such resonant 

transformation effects have also been studied in the 
context of photon-axion oscillations ■ 

However, the case of photons with different polariza- 
tion has been ignored until recently when Broderick and 
Blandford [6[ for the first time considered the Faraday 
effect across regions of slow i?-field reversal where the 
dispersion relations of photons with opposite helicities 
cross over. If this transition is adiabatic 2 , photons flip 
their helicity. The practical astronomical consequence is 



that the usual rotation measure for photon polarization 
RM ex J B • d£ changes to ± J |B • d£| and thus can build 
up continuously across regions of opposite field direction. 
The adiabaticity condition depends on photon frequency. 
Comparing the Faraday effect at frequencies above and 
below a certain u; cr ; t may allow one to study the geometry 
of astrophysical magnetic fields. 

The main purpose of our paper is to formulate the 
adiabatic Faraday effect in the more familiar language 
of particle oscillations using a linearized wave equation. 
In this way the evolution of photon polarization is 
described by a Schrodinger equation for the two-level 
system consisting of the two photon helicity states. While 
the underlying physics, of course, is precisely as discussed 
in Ref. [6], our approach helps to bring out the explicit 
analogy to neutrino oscillations and therefore allows one 
to borrow both intuition and results from more familiar 
cases. In addition, we recover the nature of Faraday 
rotation at intermediate adiabaticity. 

The article is organized as follows. In Sec.[n]we derive 
a two-level Schrodinger equation for photon polarizations 
and identify the condition for adiabatic helicity flipping. 
In Sec. IIIII we work out the impact on linearly polarized 
states and recover the adiabatic Faraday effect. We also 
discuss a signature for the transition from adiabatic to 
non-adiabatic evolution. We conclude in Sec. IIVI 



1 Optical birefringence can be used to explain neutrino oscillations 
in a pedagogical demonstration experiment [ij • 

2 Broderick and Blandford use the term "super adiabatic" in this 
context, but we prefer to keep the usual terminology. 



II. HAMILTONIAN APPROACH 

A. Linearizing the wave equation 

The evolution of photon polarization, neutrino flavor 
oscillations, and similar phenomena derive from the 
underlying wave equation for the relevant fields. Since 
dispersion plays a central role, it is easiest to assume 
harmonic time variation e -1 "' for all fields, i.e. to study 
the spatial variation of a monochromatic wave. Using 
natural units with Ti = c = 1, one finds a stationary 
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Klein-Gordon equation of the form 3 



V A 



n) a. 



(i) 



where A is a "spinor" of amplitudes describing the multi- 
component wave phenomenon. For photon propagation, 
this is the Jones vector formed by the transverse 
components of the electric field (Appendix A), whereas 
for neutrinos the wavefunction in flavor space comes into 
play. For photons, l~l is the "polarization tensor," whereas 
for neutrinos, l~l = M 2 with M being the mass matrix that 
is non-diagonal in the interaction basis. 

The problem simplifies further if the waves are 



relativistic, 



> 



where m is a typical 



eigenvalue of M 2 or of 11. In the language of refractive 
indices, the dispersion relation is written as k — nuj and 
the relativistic assumption amounts to \n — 1| <C 1. The 
propagating waves now involve a short length scale (the 
wave length) and a long one corresponding to the phase 
difference (the oscillation length). Assuming propagation 
in the z-direction, one can easily separate the fast and 
slow variation by the nominal substitution Q 

- (<9 2 +w 2 ) = (i<9 2 +cj)(i<9 2 -w) -> -2w(id z +w) , (2) 

where we have used — id z — > k « oj in one of the terms. 
This approximation replaces the Klein- Gordon equation 
with a linear equation 



medium. In the context of neutrino physics such effects 
are always vanishingly small, except perhaps near a 
supernova core, whereas for photons near propagation 
threshold, polarization-dependent reflection effects may 
play a nontrivial role. 

Likewise, the wave may suffer deflection caused by 
density variations in the transverse direction. Therefore, 
the requirement that all polarization components follow 
the same trajectory with sufficient precision poses non- 
trivial constraints [6j. Typically, relativistic propagation 
will be required because near propagation threshold 
the refractive indices between different components 
are largest and the differential deflection of different 
polarization states would be largest. 

The connection between a classical two-level equation 
and a true quantum equation as well as the precise role of 
the relativistic approximation will be explored elsewhere. 
Here we simply take advantage of the formal equivalence 
of our problem with a two-level quantum system. 



B. Photon dispersion in magnetized plasma 

Photon dispersion in a cold collisionless plasma with 
electron density n e and weak external magnetic field B 
is determined by the plasma frequency and cyclotron 
frequency that are respectively 4 



(3) 



Moreover, we are not interested in the overall phase of 
the wave, but only in phase differences between different 
components. Therefore, on the r.h.s. we may drop terms 
that are proportional to the unit matrix and find 



i<9 z A = HA. 



(4) 



The "Hamiltonian" is H = H/2uj = M 2 /2w. 

Sometimes Eq. (U) is written in terms of a parameter 
"time," playing the role of an affine parameter along the 
trajectory. The linear wave equation then manifestly 
resembles a Schrodinger equation for a two-level system 
such as a spin evolving in a magnetic field. We stress, 
however, that Eq. Q is a classical equation. Using 
natural units somewhat obscures that h does not appear 
anyway, at least not when considering the variation of 
photon polarization along a trajectory. 

The analogy to a Schrodinger equation reveals that 
the linear wave equation provides for unitary evolution 
along the beam, i.e., any photon that disappears from one 
polarization state must appear in another. Linearizing 
the wave equation has removed such physical effects as 
reflection by inhomogeneities or discontinuities of the 



3 We use sans-serif letters to denote matrices in polarization space 
and bold faced sans-serif letters for "spinors" in this space. 




3.7 x 10" 



"eV^) 172 , (5) 
Vcm J / 



=e — = 1.16 x 1(T 8 eV f — — ) , (6) 
m e \ Gauss ) 

where the elementary charge e was taken to be positive. 
The analogous contributions from ions are much smaller 
and will be neglected. 

We decompose B into a longitudinal component B» 
along the direction of propagation (^-direction) and a 
transverse one B± . The impact of the magnetic field is 
determined by the dimensionless parameters 



b\\ i = 



B 



1,-L 



lj B ujm e 
Propagating modes exist only if (Appendix A) 



(7) 



(8) 



4 In the particle-physics literature, rationalized units with 
a = e 2 /4?r ~ 1/137 are almost always employed, whereas 
in the context of plasma physics and photon propagation, 
unrationalized units corresponding to a = e 2 ~ 1/137 are used, 
assuming h = c = 1 in both cases. We follow the particle- 
physics tradition and note that a magnetic field of 1 Gauss then 
corresponds to 1.95 X 10 — 2 eV 2 . The critical field strength, 



defined by oj c 
(0.511 MeV) 2 A 



is then correctly found as B c 



mile 



■ X 10 11 eV^ = 4.4 X 10" Gauss 
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The assumption that we are dealing with relativistic 
waves thus implies that |&ii xl ^ 1- 

The Hamiltonian for the evolution of photon polariza- 
tion is found to be (Appendix A) 



2lu 



°ll 

e i2 V b 2 



-i2<p b 2 



(9) 



where we have dropped terms proportional to the unit 
matrix. We have used the helicity basis, i.e. the 
components of the spinor A are the amplitudes of the 
two circular polarization components of the wave. The 
angle p describes the direction of Bi relative to a fixed 
direction transverse to the photon trajectory. 

To linear order in B, the matrix H is always diagonal 
and even if B changes or reverses along the trajectory, 
nothing new happens. Photons that begin linearly 
polarized stay that way, except that their plane of 
polarization rotates by an angle proportional to J B\\dz. 
Therefore, we had to include the quadratic terms, 
responsible for the Cotton- Mouton effect. It is these 
terms that can flip the photon helicity. 

The same conclusion is reached by observing that a 
longitudinal i?-field is symmetric relative to rotations 
around the z-direction. Therefore, angular momentum 
along z is conserved and a longitudinal £?-field can not 
induce transitions between photons of opposite helicity. 
This can be achieved only by the transverse field. By 
the same token, transverse fields are required to achieve 
mixing with spin-0 particles such neutral pions or axions 
or with spin-2 particles such as gravitons. 



C. Adiabatic helicity flip 



perturbation provided by the transverse field. It prevents 
an actual crossing of the eigenvalues and leads to an 
adiabatic helicity reversal in the spirit of the MSW effect. 

Within the Landau-Zener approximation 0, Q, the 
probability for the state to jump over the avoided level 
crossing while preserving its helicity is 



(10) 



The adiabaticity parameter 7 compares the rate-of- 
change of the energy splitting with the oscillation 
frequency. When the unperturbed levels would cross 
(here the Bn reversal), one finds the usual result 



4|Hi 



l H 22 ~ H lll 



where a prime denotes d/dz. In our case this is 



1=7T- 



2lo \2b' 



Alo 4 



(11) 



(12) 



where we have used that at a B\\ reversal B = B± and 
therefore b± = u> c /u>. We have also introduced the length 
scale of B- variation £g — \BUB_l\- If the magnitude B 
is fixed so that B\\ = Bcos8 and B± = BsinO with 
9 the angle between photon direction and B, we find 
Iq = l^'l -1 . Moreover, using cycle frequencies v = gj/2it, 
we find that the upper v limit given in Eq. (6) of Ref. @ 
corresponds in our treatment to 7 = 1. In other words, 
7=1 defines a critical frequency 



1/4 



(13) 



If photons propagate in the presence of a purely 
transverse -B-field, the Cotton-Mouton effect provides a 
refractive difference between the modes that are linearly 
polarized parallel and orthogonal to the field. A photon 
that is initially circularly polarized will acquire elliptic 
polarization, eventually reverse its helicity, and later 
return to the initial polarization in the spirit of an 
oscillation phenomenon. 

However, a realistic situation is different. Astrophys- 
ical magnetic fields are not homogeneous and vary in 
magnitude and direction along any given photon path. In 
the relativistic limit where |&||,x| "C 1, the refractive effect 
of S11 is much larger than that of B± because the Faraday 
effect is linear in B, the Cotton-Mouton effect quadratic. 
Therefore, almost everywhere the "mass difference" of 
photons with opposite helicity is large compared with the 
mixing energy provided by the transverse field. In other 
words, the unitary transformation between helicity and 
propagation eigenstates involves a small mixing angle. 

However, Bi\ may vanish somewhere and may reverse 
while B± remains at a typical value. At the Bu inversion 
point the helicity eigenstates become degenerate and 
their effective masses cross over, except for the small 



such that we are in the adiabatic regime for uj <C w C rit 
and in the nonadiabatic regime for uj 3> uJmt- 

The Landau-Zener approximation is strictly applicable 
only when B\\ decreases linearly and certain other 
conditions are met [9J. There is however substantial 
literature on its refinements [Toj . In particular, the 
"double exponential" ansatz parameterizes the jump 
probability in a more widely applicable way [ill ]. We 
will continue to use the Landau-Zener formula for its 
simplicity, but our discussion proceeds unaltered if one 
substitutes a more accurate expression. 

Thus far we assumed that the transverse i?~field points 
in a fixed direction, but of course this direction changes 
on the same length scale lg as all other properties. 
This may have interesting consequences fl2j . However, 
in the adiabatic approximation, i?-field twisting is slow 
compared with the oscillation length and thus would not 
produce new effects. Rapid twisting would lead to non- 
adiabatic transitions. 

The B-field may have other variations, for example 
fast variations caused by turbulence. The evolution 
of photon polarization in such circumstances could be 
performed on the level of our Schrodinger equation. Both 
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for solar and supernova neutrino oscillations it was found 
that relatively small density fluctuations of the medium 
can severely affect the MSW effect [H El- Similar 
phenomena for radio waves are expected and would be 
an interesting subject of study. 



III. FARADAY ROTATION 

A. Ordinary Faraday effect 

We now investigate the role of adiabatic helicity flips 
on the rotation of linearly polarized radiation traversing 
a region of magnetized plasma. If we begin with Ao at 
the source, at the detector we will have 

A D = UtotAo . (14) 

The unitary matrix taking the initial to the final state is 



Utot = S exp -i 



Hdz 



(15) 



where the Hamiltonian in the helicity basis is given by 
Eq. ©. The parameters w p , b\\ t ± and ip all depend on 
z. At different locations the matrices H(zi) and Hfa) do 
not in general commute, so the exponential is understood 
in the space-ordering convention symbolized by S. Note 
that we work in the helicity basis. However, if photons 
are produced and detected in regions that do not involve 
large transverse fields, the helicity states are identical 
to the propagation states at source and detector. Thus 
oscillatory terms depending on the mixing angle at source 
and detector, as in Ref. [II [3, will vanish. 

In the extreme case when the transverse field B±_ is 
either zero or so small that it plays no role, H is always 
diagonal and one finds explicitly 



Ut 





3 +i< 



with the phase 



I 

2cj 



dz Wp by 



2w 2 m 2 



,B-d£. 



(16) 



(17) 



The two helicity states pick up equal but opposite 
phases, implying that the plane of polarization rotates 
by the angle <f>: the usual Faraday effect. However, the 
polarization direction at the source is not known, so the 
measurable quantity is the variation of <j> with frequency. 
This "rotation measure" is 



RM 



A 2 87r 2 m? 



n e B ■ d£. 



(18) 



where ui = 2tt/X has been used, with A being the 
wavelength. RM itself does not depend on frequency. 



B. Adiabatic Faraday effect 

In order to understand the impact of an adiabatic 
helicity flip we subdivide Utot into several pieces. 
To be specific we assume a single B-field reversal. 
The propagation up to somewhat before this point is 
described by the ordinary Faraday effect. The same is 
true after the reversal onwards. Therefore, the overall 
effect is 



with 



Ui, 



Utot = U 2 U fli pUi 



e + i ^ 1 - 2 



(19) 



(20) 



The phases are given by ordinary Faraday integrals on 
path 1, leading from the source to the reversal, and 
path 2, leading from the reversal to the detector. If the 
reversal is non-adiabatic, Ua; p is the unit matrix and we 
recover the previous result: 4> — 4>i + (f>2- 

On the other hand, if the reversal is perfectly adiabatic 
it has the effect of exchanging the left- and right-handed 
helicity states 



Uflip 



1 

1 



(21) 



In general, an additional relative phase arises — we will 
be more precise later. The flip matrix implies 



Utot = 



e K0i-02) 

g-i(01-02) 



(22) 



If the initial linear polarization state is Ao = (1, 1) the 
non-adiabatic (na) and adiabatic (ad) final states are 



«na 



1 + i(0 1 +0 2 ) 



i ad 



e -KOi-02) 

e -i(0l-0 2 ) 



(23) 



The impact of the second part of the trajectory is the 
same in both cases: the helicity components acquire a 
relative phase <f>2, whereas the impact of the first part is 
reversed. The overall rotation of linear polarization is 



2 and $ ad = -0i + 02 ■ 



(24) 



In other words, the adiabatic helicity flip has the effect of 
reversing the rotation measure accrued on the path before 
the flip. This is physically understood: The helicity 
states in terms of electric fields are [E x ± iE y )/^/2. 
Exchanging helicities amounts to E y — > —E y and thus 
reverses the polarization angle relative to the x-direction. 

We juxtapose the non-adiabatic and adiabatic evolu- 
tion of RM in Fig.[T]that is analogous to Fig. 2 of Ref. [6j. 
An important difference is that for the adiabatic case 
these authors actually show the evolution of the phase 
difference between the two propagation eigenmodes, a 
quantity that indeed always increases. However, the only 
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non-adiabatic 



adiabatic 



FIG. 1: Evolution of the rotation measure in a magnetized 
plasma with changing B-field orientation. The ordinary (non- 
adiabatic) and adiabatic cases are shown. 



observable quantity is the RM. It acquires a minus sign 
at each adiabatic reversal and thus jumps by a large 
amount. The overall adiabatic RM is 



RM 



ad 



± 



87r 2 m^ 



/ n e \B-d£\ 
Jo 



(25) 



The absolute sign is identical with the sign of Bn on the 
last sub-trajectory. 



At the fully adiabatic or fully non-adiabatic limit, 
sin 2£ = and the two helicities arrive with equal 
intensity: the beam is linearly polarized. At intermediate 
adiabaticity this is no longer true and the final beam is 
elliptically polarized to a degree depending on frequency. 
Using Eqs. (JT2J) and (jl3]l . we find 



cos £ = cxp 



-)' 

^crit / 



(29) 



whereas from Eq. (|16[) we find that (f>i varies as 
J dzujp &||/2cj, the integral extending from the source to 
the point of reversal. 

The degree of elliptical polarization thus has two 
variations imprinted upon it as a function of frequency. 
One arises from the usual RM accrued between source 
and field reversal which determines the orientation with 
which the linear polarization enters the cross-over region. 
The other is a slow variation determined by the jump 
probability as a function of frequency, approximately 
given by the Landau-Zener formula. If one could 
measure the polarization state over a reasonably broad 
range of frequencies, a measure of -B-field tomography 
would become possible. Unlike the adiabatic or non- 
adiabatic limit, at intermediate adiabaticity one has 
the opportunity to probe the magnetic fields in specific 
segments of the path traveled. 



C. Intermediate adiabaticity 

The final polarization of the beam depends crucially 
on the degree of adiabaticity. The entire effect of the 
reversal is parameterized in the flip matrix. In general, 
Uflip is given by a unitary matrix 



U 



flip 



cos£ 



3 i<5 2 



sin£ 



sin£ e~ ldl cos£ 



(26) 



where cos£ = yJ~P- } and Pj is the probability to jump 
between helicities at the crossing. Note the presence of 
other phases acquired at the reversal. The non-adiabatic 
limit is obtained for cos£ = 1 with other phases set 
to zero. Similarly the adiabatic limit is obtained for 
sin£ = 1, which is the same as in Eq. ([2Tj) for overall 
phase choices a = 82 = tt/2. 

If the initial linear polarization state is Ao = (1,1), the 
final state (up to an overall phase) is 



cos£ 



a i(<52+20i-02) 



sin£ 



-i(5i-20i-0 2 ) 



cos£ 



-i(<5 2 — <fe) 



sin£ 



(27) 



The intensity of the left and right helicity components of 
this state are 



I± = l±sin2£cos((5-20i) , 



(28) 



where 8 — 8\ — 82 is a combination of the phases acquired 
at the reversal. 



D. Multiple reversals 

We have mostly confined our discussion to a single field 
reversal, but the formalism is easily extended to multiple 
reversals. As emphasized by Broderick and Blandford, 
the ordinary Faraday effect for N 3> I domains, each 
causing a rotation fa, adds up to 



N 



cx VN 



(30) 



On the other hand, the adiabatic Faraday effect leads to 
a rotation by 



N 

$ ad = ^ fa cx N 



(31) 



Clearly, for media with multiple reversals this leads to a 
larger Faraday rotation at adiabatic frequencies. 

For quasi-adiabatic frequencies, multiple reversals lead 
to loss of linear polarization by each of the reversals. 
Schematically, for N domains with N — 1 field reversals, 
the intensities of the helicity components arriving at the 
detector is 



1+ 
I- 



N—l 

n 



cos 2 & 
sin 2 & 



sin 2 f j 
cos 2 & 



(32) 
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where we have ignored the fast oscillatory terms. 
Clearly, there is increased elliptic polarization at the 
critical frequencies corresponding to each of the magnetic 
field reversal. A more precise prediction, including 
interference effects due to the oscillatory terms, can be 
calculated using the recipe prescribed in Ref. fl7j ]. 



IV. CONCLUSIONS 

We have studied the adiabatic Faraday effect that was 
recently discovered by Broderick and Blandford. We have 
used a simple formalism that linearizes the photon Klein- 
Gordon equation and amounts to a Schrodinger equation 
for a two-level system consisting of the photon helicities. 
This approach is commonplace in the context of neutrino 
flavor oscillations. Once the photon dispersion relation 
has been identified, the adiabaticity condition follows 
immediately from well-known textbook results such as 
the Landau-Zener approximation. Formulating the 
adiabatic Faraday effect in a language familiar from 
neutrino physics may allow for a broader appreciation 
of this intriguing "MSW effect for photons." 

Despite the similarities, Faraday rotation also has 
important differences with neutrino oscillations. For 
neutrinos, the initial states are always weak-interaction 
eigenstates, whereas photons can in general be produced 
in any polarization state. During propagation, neutrinos 
typically have non-maximal mixing except on the MSW 
resonance, causing large flavor transitions. For photons, 
the plane of polarization is typically rotated by many full 
cycles independently of any transverse field. Therefore, 
the adiabatic helicity flip does not enhance the transition 
between linear polarization states, it modifies the 
rotation measure, the way the rotation of the plane 
of polarization varies with frequency. Finally, unlike a 
neutrino detector, a photon detector in the radio band 
has the capability of identifying not merely the intensities 
of each polarization state, but also the relative phase 
between them. So the observational manifestation of the 
MSW effect is rather different for neutrinos and photons. 

To summarize, the adiabatic Faraday effect is the 
process of photon helicity getting flipped adiabatically 
in regions of slow B-field reversal, and has the effect of 
reversing the rotation measure accrued by the photon up 
to that point on its trajectory. Therefore, the subsequent 
Faraday effect goes effectively in the same direction as the 
Faraday effect before the reversal region, giving much 
larger overall rotation measure. The global sign of the 
rotation measure depends only on the i?-field direction 
on the last leg of the photon path. 

The rotation measure both in the adiabatic and 
non-adiabatic regime is independent of frequency. 
Observations over a broad range of frequencies, ranging 
from below to above the critical frequency, would 
reveal a transition range of fast modulations of elliptic 
polarization. In the adiabatic regime, similar variations 
can arise in the transverse direction across the sky. 




non-adiabatic 



cut-off 



(0. 



CO 



crit 



FIG. 2: The rotation measure as a function of frequency. The 
ordinary (non-adiabatic) and adiabatic regimes are shown. 



The astronomical potential of these effects for B- 
field tomography has been explored by Broderick and 
Blandford. The main purpose of our note was to clarify 
the basic principles of the adiabatic Faraday effect that 
is a neat application of the formalism usually applied to 
the MSW effect in neutrino oscillations. 



Note Added in Proof 

After our manuscript had gone to press, D. Melrose 
circulated a preprint [fij ]. claiming that the adiabatic 
modification of the Faraday effect was not observable. 
We believe that underlying this critique is a misunder- 
standing. 

We fully agree with the formalism used by Melrose. 
Our complex two-spinor A is equivalent to a 2 x 2 
density matrix Sy = A^A* and our Eq. ([?]) can be 
written as a commutator equation in the form i<9 z S = 
[H,S]. Moreover, any Hermitean 2x2 matrix R can 
be written in terms of a three-vector R in the form 
^Tr(R) + iR ■ er with er being a vector of Pauli matrices. 
Our Eq. (Q| is equivalent to d z S = H x S where S is 
the vector representing S and H the one representing 
H. Therefore, our Eq. (j4|) is equivalent to the "spin- 
precession equation" described after Eq. (4) of Ref. [19| . 
In particular, we agree, that S 2 is invariant and that 
H • S is an adiabatic invariant if H changes slowly as 
a function of z. A precessing spin following a slowly 
changing B-field is the usual visualization of adiabatic 
neutrino oscillations. 

The adiabatic limit is defined by the oscillation 
length being shorter than the length scale of field 
reversal, implying that the position angle (PA) of linearly 
polarized light must undergo many Faraday revolutions 
between source and detector. Therefore, we agree 
that PA, being measurable only modulo 7r, carries no 
information about the difference between normal and 
adiabatic Faraday effect. This information is contained in 
how quickly PA changes with frequency, i.e. the rotation 
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measure (RM). In our Fig. Q] we have sketched how RM 
varies as a function of distance of a hypothetical observer 
from the source. However, a real observer is at a fixed 
location and observes the effect as a function of frequency 
as sketched in Fig. [2] The observed RM in the normal 
and adiabatic regime can be different, and are given 
by Eqs. (TT5|) and (l2~5j) respectively. Of course, in the 
transition region between adiabatic and ordinary Faraday 
effect, the photon polarization is elliptical. Eq. (|28|) 
shows that unless the probability of helicity flip is zero 
or one, the light is always elliptically polarized . Since 
the RM is not strictly defined in that case, we sketch the 
RM for the major axis of the polarization ellipse. 
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The lowest frequencies arise in the homogeneous case 
when all spatial derivatives vanish and the l.h.s. of 
Eq. (|A1|) vanishes identically, leading to 



u> = Up or u> = | lo c + y \u>l + u>l . (A3) 
The first solution corresponds to the ordinary plasma 
oscillation. The other connects to the dispersion relation 
for propagating modes for non-vanishing k and thus is 
the minimum frequency required for a propagating wave. 

For propagating modes the ^-equation of Eq. (|A1[) 
represents a constraint, allowing one to eliminate E z . 
We represent the propagating modes by the usual Jones 
vector A = (E x ,E y ), leading to a stationary Klein- 
Gordon equation in the form 



d 2 A = (u 



uj 2 -cu'*n ) a 



(A4) 



The 2x2 reduced polarization matrix is found to be 



n« = ii,.. 



n tz n ZJ 



, i,j = x or y. (A5) 



With the dimensionless parameter b p — uj p /uj we find for 
n the expression 



Appendix A: Photon dispersion in a cold 
magnetized plasma 

The propagation of electromagnetic waves in plasma is 
governed by Maxwell's equations, including polarizations 
and currents induced by the response to the wave. 
Assuming a homogeneous cold magnetized plasma, 
harmonic time variation e _ILJt for all quantities, and 
spatial variation of the wave only along the z-direction, 
one finds a stationary Klein-Gordon equation for the 
electric field vector E of the form [T|| 



'ci 2 



E = Uj 



E. 



(Al) 



where uj p is the plasma frequency. II is the reduced 
polarization tensor (polarization tensor in units of uj 2 ) 
which depends on uj and the medium properties. 



Assuming the B field has components B z 



-B|| , B v 



B± and B y = and using the b parameters defined in 
Eq. one finds 



n 



i 



l-b 2 



b 2 ± 



-i&U 
-hb i 



i&H 
f 

-ib. 



-&H&L \ 



(A2) 



f 



l-b 2 p -b 2 ± i&||(l-6 2 ) N 



(l_&3)(l-6jj)_&a V-^l-fcp) l - h \ 



(A6) 

Notice that the absence of absorption renders this matrix 
Hermitean. 

If all of &p,j_,|| are small compared to unity and we 
expand up to quadratic order we find 




In general the B-field has components B± cos <p, B± sin ip 
and -By relative to a coordinate system with a fixed re- 
direction, implying that the polarization matrix must be 
rotated correspondingly. Moreover, for some purposes it 
is simpler to work in the helicity basis where the electric 
field has the two components (E x ±iE y ) / ' Altogether 
these transformations lead to 



fl = 1 + b\ + \b\ 



°ll 

■ie i2 ^ b\ 



2 e o_,_ 



(A8) 



This is the result used in the main text. 
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